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£SJ , Abstract. The largest uncertainty in calculations of hydrogen ground-state hyperfine 

splitting comes from corrections due to proton stucture. We review these corrections, 
^ ' with special mention of the inelastic, or polarizability, corrections which have been 

recently re-evaluated. Summing up the arguably best current values for the calculated 
corrections leaves us 1-2 ppm short of the experimental data. We speculate how one 
I/") ' may explain this shortfall, and along the way broadly outline the derivations of the 

relevant formulas, attempting to explain how certain terms come to appear and com- 
menting on the use of unsubtracted dispersion relations. 
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1 Introduction 

Hyperfine splitting (hfs)in the hydrogen ground state is measured to 13 signifi- 
cant figures in frequency units pQ, 



E his (e-p) = 1 420.405 751 766 7(9) MHz . (1) 

Oh. Theory is far from this level of accuracy, and theorists are hopeful of obtaining 

O 1 calculations accurate to a part per million (ppm) or so. We are close to reaching 

D , this goal, but some improvement is still needed and there currently seems to be 

a few ppm discrepancy between the best calculations and the data. 

The main uncertainty in calculating the hfs in hydrogen comes from the 
hadronic, or proton structure, corrections. One can contrast this to the case 
of muonium, where the "nucleus" is a point particle, so that calculations are 
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d ' almost purely QED, and agreement between theory and experiment is about 0.1 

ppm [I]. 

For ordinary hydrogen, as we have said, one must consider the proton struc- 
ture, and find that it contributes about 40 ppm to the hfs. Working out these 
contributions theoretically requires knowing details about proton structure that 
cannot currently be obtained from ab initio calculation. Instead, one has to 
measure information about proton structure in other experiments, particularly 
experiments on elastic and inelastic electron-proton scattering. Then calcula- 
tions are done to relate the scattering information to the bound state energy. 
Recent new results have been driven by improvement in the data, including both 
new data for polarized inelastic scattering in kinematic regions of interest to hfs 
calculators and new analyses of the elastic scattering data. These will be outlined 
below. 
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Historically, the elastic and inelastic contributions, the latter also called po- 
larizability corrections, have often been treated separately, with the elastic cor- 
rections further divided into a nonrelativistic Zemach term and relativistic re- 
coil corrections. From a modern viewpoint, the elastic and inelastic corrections 
should be treated as a unit since the sum lacks certain ambiguities that exist 
in the individual terms. The present discussion will focus on the polarizability 
contributions, but following the last remark, discussion of the Zemach and recoil 
corrections will not be omitted. 



2 Hyperfine splitting calculations 

The calculated hyperfine splitting in hydrogen is [l|2|3j 

E hfs (e-p) = (1 + Zi QED + Z\^ cak + A StI ) E p F , (2) 
where the Fermi energy is 

8a 3 ml 8a 4 m 3 

with m r — m e m P l '(m p + m e ) being the reduced mass (and there are hadronic 
and muonic vacuum polarization terms [2] which are included as higher order 
corrections to the Zemach term below). The QED terms are accurately calculated 
and well known. They will not be discussed, except to mention that could be 
obtained without calculation. The QED corrections are the same as for muonium, 
so it is possible to obtain them to an accuracy more than adequate for the present 
purpose using muonium hfs data and a judicious subtraction [415) . The weak 
interaction corrections [6] also will not be discussed, and are in any case quite 
small. We will discuss the proton structure dependent corrections, 

zistr - A el + A inel = A Z + A P R + A pol , (4) 

where the terms on the right-hand-side are the Zemach, recoil, and polarizability 
corrections. 

Generically, the proton structure corrections come from two-photon exchange, 
as diagramed in Fig. [1] The diagram can be seen as Compton scattering of off- 
shell photons from an electron knit together with similar Compton scattering 
from the proton. (We are neglecting the characteristic momentum of the bound 
electron. This allows a noticeably simpler two-photon calculation than for a 
scattering process [7 . One can show that keeping the characteristic momentum 
would give corrections of 0(am e /m p ) smaller than terms that are kept [5].) 



2.1 Elastic terms: the Zemach correction 

In this author's opinion, the best calculation of the box diagram uses dispersion 
relations, even though there are questions about using dispersion relations in 
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Fig. 1. Generic two-photon exchange diagrams, giving proton-structure corrections to 
hyperfine splitting. 



their unsubtractcd form. In a dispersive calculation, it is easy to consider the 
elastic and inelastic intermediate states simultaneously. This one should do be- 
cause the full calculation is well defined, even though historically terms have 
been shuttled between the "elastic" and "inelastic" contributions. 

At the outset, however, we will present results from a direct calculation of the 
elastic contributions, without dispersion theory. The results have been obtained 
by a number of authors [9110] , and follow after assuming a certain photon-proton- 
proton vertex which is plausible but which cannot be defended perfectly. 

The "elastic' contibutions are those where the hadronic intermediate state, 
the blob in Fig. [TJ is just a proton. The diagram specializes to Fig. The 
photon-electron vertex is known, and we use [9110] 

r, = 7„*i (i 2 ) + ^—o^F^q 2 ) (5) 

2m p 

for the photon-proton vertex with incoming photon momentum q. Functions F± 
and F2 are the Dirac and Pauli form factors of the proton, which are measured in 
elastic electron-proton scattering. The normalization is i*i(0) = 1 and ^(0) = 
K p , where n p is the proton's anomalous magnetic moment measured in proton 
magnetons. 

The above photon-proton vertex is complete and correct only if the protons 
entering and exiting the vertex have physical, on-shell, momenta. In a loop di- 
agram, the intermediate proton is on-shell only at special values of momenta 
out of a continuum of momenta. Hence one may feel some hesitation in using 
the results that follow. However, there is also the (coming) dispersive calcula- 
tion which only needs the vertices when all protons are on-shell, and so (again 
modulo questions surrounding dispersion relations with no subtraction) gives a 
reliable result. Any terms in the "elastic" calculations that appear to require 
modification can be fixed by adding or subtracting terms in other parts of the 
quoted result. 

The elastic contributions are separated as 

= Az + A p R . (6) 

el 

The separation is into non-relativistic and relativistic terms — "Zemach" and 
"recoil." Non-relativistic means the limit m p — ► 00 with m e held fixed and with 
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Fig. 2. Two- photon exchange diagrams for the elastic proton-structure corrections to 
hyperfine splitting. 



proton size held fixed; proton size information is embedded in the form factors 
Fi and F 2 . 

The Zemach correction was worked out by Zemach in 1956 in modern 
form it is 



8am r f 00 dQ 



2 



GE( _ Q . ) GMt&)_ 1 

t + Kp 



-2am r rz , (7) 



where the last equality defines the Zemach radius rz and we have used Q 2 = —q 2 . 
The charge and magnetic form factors are linear combinations of F\ and F%, 



Gm — Fi + F: 



2 ■ 



Q 2 

Ge = Fi — - — ^F 2 . (8) 

Table [1] gives the evaluated Zemach radius rz and correction Az for two 
modern form factor fits, and for the dipole fit, which is out of date and included 
only because it is a common benchmark. We believe the Ingo Sick fit is best for 
the purpose at hand, because the Zemach integrals depend mainly on the form 
factors at low Q 2 , and Sick's fit concentrates on the low Q 2 scattering data. 



Table 1. Values of the Zemach radius and the Zemach corrections for selected fits 
to proton elastic form factors. (The Zemach term Az includes a 1.53% correction 
from higher order electronic contributions |12j . as well as a +0.07 ppm correction 
from muonic vacuum polarization and a +0.01 ppm correction from hadronic vacuum 
polarization [2].) 

Form factor rz (fm) Az (ppm) 

Kelly [I3] 1.069(13) -40.93(49) 

Sick [13] 1.086(12) -41.59(46) 

dipole 1.025 -39.24 



2.2 Elastic terms: recoil corrections 

The relativistic elastic corrections A P R are known as recoil corrections. They 
depend on the form factors and hence are part of the proton structure corrections. 
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However, evaluating A P R with different form factor representations based on fits 
to the scattering data reveals that its numerical value is fairly stable (to about 
±0.15 ppm) by present standards. 
The full result is 

_ am e m p 

R ~ 2(1 + K p )n(ml - ml) 



dQ 2 fl3i{T p )-A^r- f3 1 (r e )-4 y 
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Q 2 

dQ 2 



Fi(-Q 2 )Ga/(-Q 2 ) 



(/3 2 (t p ) - /3 2 (r e )) F 2 (-g 2 )G M (-Q 2 ) 



lo Q 

2am r f°° „2,Gm{-Q 2 ) 



irm 2 



Jo t + K p 



am e r°° dQ 2 

2(1 + K p )nm p J Q 2 



^{T e )F 2 {-Q 2 ) , (9) 



where /3i t 2 are auxiliary functions that were first found useful in discussing the 
inelastic terms |14|15|16|17plp9] . 



1 (t) = -3r + 2r 2 + 2(2 - t)^t{t + 1) = -/3(t) , 



/3 2 (r) = l + 2r-2vV(r+T). (10) 
These are used with the notation, 

(11) 



4r7ip ' 4m 2 

The reasons for showing the whole formula for the recoil corrections is partly 
to show it is not so long (it sometimes appears more forbidding, Cf. [9]), to 
explicitly see the form factor dependence, and to display the F 2 term in the last 
line. The F% term is noteworthy because it is absent in a dispersive calculation, 
in contrast to all the other terms, which come forth unchanged. 

All the integrals are finite, although some require the form factors to prevent 
ultraviolet divergence, and some would be infrared divergent if m e — * 0. Of minor 
interest, the penultimate line could be subsumed into the Zemach correction, if 
the Zemach correction were to be written in terms of F%Gm rather than GeGm- 



2.3 Inelastic terms: polarizability corrections 

When the blob in Fig. [TJ is not a lone proton, we obtain inelastic contributions 
or polarizability contributions [14115116117118119] . The inelastic contributions are 
not calculable ab initio. Instead, one relates them to the amplitude for forward 
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Compton scattering of off-shell photons off protons, given in terms of the matrix 
element 

T^(q,p, S) = ^— J e*** ( P S\ j M (£)j,(0) \pS) , (12) 

where is the electromagnetic current and the states are proton states of mo- 
mentum p and spin 4-vector S. The spin dependence is in the antisymmetric 
part 



rpA " O 



I 



(13) 



p-q 

There are two structure functions Hi and H 2 which depend on q 2 and on the 
photon energy v, defined in the lab frame so that m p v = p ■ q. 

There is an optical theorem that relates the imaginary part of the forward 
Compton amplitude to the cross section for inelastic scattering of off-shell pho- 
tons from protons. The relations precisely are 

[ jTn 

ImH^q 2 ) = - gi{v,q 2 ) and Im H 2 (u, q 2 ) = -£ g 2 {v, q 2 ) , (14) 

where gi and g 2 are functions appearing in the cross section and are mea- 
sured [20121122123124] at SLAC, HERMES, JLab, and elsewhere. 

Using the Compton amplitude in terms of Hi and H 2 , Eq. (fl~3|). in evaluating 
the inelastic part of the two-photon loops gives 

2am e f d A Q 



^pol — "FT - 



el (l + K p )TT 3 m p J (Q 4 + Am 2 e Q 2 )Q 2 
{(2Q 2 + Q 2 )H\ nel (iQ 0) -Q 2 ) - 3Q 2 Q 2 H™ el (iQ 0) -Q 2 )} , 



(15) 



where we have Wick rotated the integral so that Qo = —iu, Q = q, and Q 2 = 
Q 2 +Q 2 . Since Hi 2 are n °t measured, we obtain them from a dispersion relation, 
which will discussed in a subsequent section. Assuming no subtraction, 

rrinelf 2, 1 f° , l^Hi{ V ,q 2 ) 

H l ("1,Q ) = - L dv ~2 _ v 2 > ( 16 ) 

where the integral is only over the inelastic region (yth = m^ + lm 2 +Q 2 )/(2m p )), 
and similarly for in- 
putting things together, neglecting m e inside the integral, and integrating 
what can be integrated, one obtains the expression 

^'= 2(l + I>,,^ ' + ^ <"> 
where, with r = v 2 /Q 2 , 

Ai = lJ™^{Fi(-Q 2 ) + 4 mp J™^P(T)gi(v,-Q 2 )} , (18) 
Ao = -12 



f°° dQ 2 [°° dv 
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The integral for A\ is touchy. Only the second term comes from the procedure 
just outlined. The first arises when one applies the dispersive calculation also to 
the elastic corrections, and discovers the F 2 2 term pointed out earlier in Eq. ((9]) is 
absent. It was then thought convenient to add the first term as seen above, and 
then subtract the same term from the recoil contributions. This leaves the elastic 
corrections exactly as already shown. This stratagem also allows the electron 
mass to be taken to zero in A\. The individual dQ 2 integrals in A\ diverge (they 
would not had the electron mass been kept), but the whole is finite because of 
the Gerasimov-Drell-Hearn (GDH) sum rule |25|26j . 



coupled with the observation that the auxiliary function (3(t) becomes unity as 
we approach the real photon point. 

The polarizability expressions have some history. A short version is that 
considerations of A po \ were begun by Iddings in 1965 [14] . improved by Drell 
and Sullivan in 1966 [15j, and given in present notation by de Rafael in 1971 [16] . 
But no sufficient spin-dependent data existed, so it was several decades before the 
formula could be evaluated to a result incompatible with zero. In 2002, Faustov 
and Martynenko became the first to use g\ y 2 data to obtain results inconsistent 
with zero [IS]. Their 2002 result was 



However, they only used SLAC data and A\ and Ai are sensitive to the behav- 
ior of the structure functions at low Q 2 . Also in 2002 there appeared analytic 
expressions for fit to data by Simula, Osipenko, Ricco, and Taiuti [27] . 
which included JLab as well as SLAC data. They did not at that time inte- 
grate their results to obtain A po \. Had they done so, they would have obtained 

/V = (0.4±0.6) P pmQl. 

We now have enough information to discover a bit of trouble. Table [2] sum- 
marizes how things stood before the 2005/2006 re-evalations of A po \. The sum of 
all corrections is 1.59±0.77 ppm short of what would be desired by experimental 
data. Using the Simula et al. value for A po \ would make the deficit greater. Using 
other proton form factor fits (limitng ourselves to modern ones that fit the data 
well) in evaluating Az can reduce the deficit somewhat, but not by enough to 
ameliorate the problem [19] . 

The discrepancy is not large, measured in standard deviations. On the other 
hand, the problem is clearly not in statistical fluctuations of the hfs measurement 
one is trying to explain, so one would like to do better. As listed in the Table, the 
largest uncertainty in the corrections comes from Z\ po i. Further, the polarizability 
corrections require knowledge of g\ and (72 at relatively low Q 2 , and good data 
pressing farther into the required kinematic regime has relatively recently become 
available from JLab (the Thomas Jefferson National Accelerator Laboratory, in 
Newport News, Virginia, USA). Accordingly, we shall present a state-of-the- 
art evaluation of the polarizability correction for electronic hydrogen. To give 




(19) 



A pol = (1.4 ±0.6) ppm 



(20) 
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Table 2. Corrections to hydrogenic hyperfine structure, as they could have been given 
in 2004. The first line with numbers gives the "target value" based on the experimental 
data and the best evaluation of the Fermi energy (8 figures) based on known physical 
constants. The corrections are listed next. (The Zemach term includes a f .53% correc- 
tion from higher order electronic contributions [12j . as well as a +0.07 ppm correction 
from muonic vacuum polarization and a +0.01 ppm correction from hadronic vacuum 
polarization [2].) The total of all corrections is 1.59+0.77 ppm short of the experimental 
value. 



Quantity 


value (ppm) uncertainty (ppm) 


(E his (e-p)/E p F ) - 1 


1 103.49 


0.01 


Aqed 


1 136.19 


0.00 


A z (using Friar & Sick gS]) 


-41.59 


0.46 




5.84 


0.15 


A pol (from Faustov & Martyenko, 2002 [18]' 


) 1.40 


0.60 


weak 


0.06 




Total 


1101.90 


0.77 


Deficit 


1.59 


0.77 



away our results [19] at the outset, we essentially confirm (remarkably, given the 
improvements in data) the 2002 results of Faustov and Martyncnko. 

3 Re-evaluation of A po \ 

Data for gi(v, q 2 ) has improved due to the EG1 experiment at JLab,, which 
had a data run in 2000-2001. Some data based on preliminary analysis became 
available in 2005 [24] ; final data is anticipated in late 2006. A sample of the new 
data is shown in Fig. [3] Since a function of two variables can be complicated to 
show, what is shown is the integral 

f* OO 7 

7(Q 2 )EE4m p / ^ 9l (v-Q 2 ) 7 (21) 

Jvth V 

which differs from an integral appearing in A\ in lacking the auxiliary function. 
The integration was done by the experimenters themselves. We remind the reader 
that this integral is expected to reach the Gerasimov-Drell-Hearn value — k 2 at 
the real photon point, and that because of cancellations the difference of this 
integral from — k 2 is more relevant to the final answer for Z\ po i than its absolute 
value. 

All the data shown is new; there was no polarized electron-proton scattering 
data available below Q 2 = 0.3 GeV 2 when Simula et al. and Faustov and Marty- 
nenko did their earlier fits. A curve obtained by integrating Simula et a/.'s fit for 
gi is also shown on the Figure; we do not have enough information to produce 
a similar curve for Faustov and Martynenko's fit. 

Integration in the region Q 2 > 0.05 GeV 2 is done using analytic fits to 
actual data for g%. For Q 2 below 0.05 GeV 2 , where there is no data, we do 
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0.05 0.1 0.15 0.2 0.25 0.3 



Q 2 (GeV 2 ) 

Fig. 3. Data for (an integral over) the spin-dependent struture function gi(v,q 2 ). The 
data is from the EG1 experiment at Jefferson Lab and from the year 2005. The Simula 
et al. curve is from a fit published in 2002. 



an interpolation based on a low Q 2 expansion within the integral to get (with 



> 0.05 GeV 2 ), 



4m, 



-~r 2 P K 2 p + 18m;ci ^ 



5mz 



7o 



%P{r) gi { V ,-Q 2 



Ql + 0(Q\) . 



(22) 



Here rp is from the expansion of the Pauli form factor F2(—Q 2 ) = k 2 (l — r p Q 2 /6 + 
and the "forward spin polarizability" 70 has been evaluated from data for other 
purposes [2"9"j . 



7o 



2a 



dv 



gi(f,0) = [-1.01 ±0.08 ±0.10] x 10" 4 fm 4 



(23) 



"P Jvth 



The parameter c\ is defined from the slope at low Q 2 of the integral shown in 

pig. m 



vth " 



I{Q 2 ) = 4m p 



8m 2 pCl Q 2 + 0(Q 4 ) 



(24) 



we find and use a = 2.95 ±0.11 GeV" 4 [19]. 

We need to comment that for A 2 , we need and there is almost no data 
for <7 2 on the proton. One estimates 172 by relating it to gi using the Wandzura- 
Wilczck relation [3D], which we shall not detail here. Fortunately, the auxiliary 
function ^(t) is small over the region where we need to do the integrals, so that 



10 Carl E. Carlson 

even when we assigned 100% error bars to the contribution from 172, the effect 
on the final answer was not great. 
Our overall result is [19] 



which is similar to the 2002 Faustov-Martynenko result. This result means that 
the polarizability corrections no longer give the largest uncertainty in Table [2 
It also means that the theory deficit outlined in Table still remains, even 
becoming modestly larger with a smaller uncertainty limit, at (1.69±0.57) ppm. 

Faustov, Gorbacheva, and Martynenko [31] quite recently published a new 
analysis and result for zi po i, obtaining the somewhat larger value 



We still believe our published result [H] is the best one for now because the 
Jefferson Lab EGlb data, which goes to lower Q 2 than other data sets, has been 
used to constrain and validate the fits that we use to do the integrals. Faustov 
et al. used only higher Q 2 data from other laboratories. It is, of course, possible 
that the final EGlb data will lead to some change. 

4 Comments on the derivations of the formulas 

The polarizability corrections depend on theoretical results that are obtained 
using unsubtracted dispersion relations. This has been alluded to before in this 
text, and this section will attempt to explain how a dispersion calculation works 
and what an unsubtracted dispersion relation is. Also, given that there may be 
a small discrepancy between calculation and data, one would like to assess the 
validity of unsubtracted dispersion relations. 

Also, the hyperfine splitting in muonic hydrogen may be measured soon. 
The polarizability corrections have been calculated for this case also [32], albeit 
only with older fits to the structure function data and the relevant formulas, 
with non-zero lepton mass everywhere, are available [32j from a single source, 
so one would like to verify these formulas. It turns out that keeping the lepton 
mass does not greatly increase calculational effort or the length of the formulas, 
so we can do the groundwork for the muonic hfs case simultaneously with the 
assessment of the ordinary hydrogen hfs calculation, although we shall not here 
display the formulas for non-zero lepton mass. 

The calculation begins by writing out the loop calculation using the known 
electron vertices and the definition of the Compton scattering amplitudes in- 
volving Hi and H2 as given in Eq. (fT3|) . One can and should use this formalism 
for all the hadronic intermediate states, including the single proton intermediate 
states. The single proton intermediate states give contributions to Hi and H2 
that can be (more-or-less) easily calculated given a photon-proton-proton vertex 
such as Eq. ([5]). For reference, we give the result for Hi, 



Ap i = 1.3 ± 0.3 ppm, 



(25) 



Z\ P oi = 2.2 ±0.8 ppm. 



(26) 



Hf = - 



2m, 




q 2 F 1 (q 2 )G M (q 2 ) , F 2 (q 2 ) 



) 



(27) 



{q 2 + ie) 2 — Am 2 v 2 Am 2 
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The criticism of the proton vertex used to obtain the above result is that it is 
not demonstrably valid when the intermediate proton is off shell, so the above 
expression may or may not be correct overall. However, it is correct at the proton 
pole. 

One may do a unified calculation of the elastic and inelastic contributions. 
Since we don't have a direct calculation of the Hi for the inelastic case, we have 
to obtain them using dispersion relations. Also obtaining the elastic terms from 
the dispersion relation is no problem |14|15j . One just needs the imaginary parts 
of H f ; these are easy to obtain, and contain Dirac delta- functions that ensure 
the elastic scattering condition v = ±Q 2 /(2m p ) and hence depends only on the 
reliable part of Eq. (|2T)) . 

Dispersion relations involve imagining one of the real variables to be a com- 
plex one an then using the Cauchy integral formula to find the functions Hi at 
a particular point in terms of an integral around the boundary of some region. 
In the present case we "disperse" in i/, treating q 2 as a constant while we do so. 
Three things are needed to make the dispersion calculation work: 

• The Cauchy formula and knowing the analytic structure of the desired am- 
plitudes. 

• The optical theorem, to relate the forward Compton Im Hi to inelastic scat- 
tering cross sections. 

• Legitimately discarding contributions from some oo contour, if the dispersion 
relation is to be "unsubtracted." 

The first two are not in question. 

For the present case, the contour of integration is illustrated in Fig. [4j where 
one should imagine the outside circle having infinite radius. The result for Hi 
begins its existence as 



The numerator of the first term is the residue (Res) from the poles in v for 
the elastic part of Hi, as from Eq. ([27]). Note that the F% term in Hf l , Eq. (|27|) . 
is constant in v, certainly lacking a pole in v. Hence this term never enters 
the dispersion relation, and no F% term arises from the elastic contribution, as 
calculated this way. 

The second term leads to the gi term in the quantity Ai given earlier, after 
using the optical theorem to relate Im. Hi to gi. 

The third term is the integral over the part of the contour which is the infinite 
radius circle. The commonly quoted results for Z\ po i, which appear in this talk, 
depend on dropping this term. The term is zero, if Hi falls to zero at infinite \v\. 
Assuming this is true, however, appears to be a dramatic assumption. It fails for 
Hi alone. Hence, for the assumption to succeed requires an exact cancelation 



Hi( V ,q 2 ) 



ResHi(v, q 2 ) 



el 



(28) 
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Fig. 4. Contour in complex v' 2 plane for applying Cauchy identity to Hi or H^- 



between elastic and inelastic contributions (or a failure of Eq. (1771) on the big 
contour). On the positive side are several considerations. One is that nearly 
the same derivation gives the GDH sum rule, which is checked experimentally 
and works, within current experimental uncertainty (8%) [33j . Also, the GDH 
sum rule is checkable theoretically in QED, where lowest order and next-to- 
lowest order perturbation theory calculations can be done, and the GDH sum 
rule appears to work 34 35 . Finally, Regge theory suggests the full Compton 
amplitude does fall to zero with energy [36], as one would like, although Regge 
theory famously gave wrong high v behavior for spin-independent analogs of gi 
and g 2 |37j . Hence there are indications, though not decisive proof, supporting 
the unsubtracted dispersion relation. 

The derivation finishes, as noted earlier, by subtracting a term involving F$ 
from the relativistic recoil term, so as to obtain exactly the elastic corrections 
A e i = Az + A P R that were obtained (say) by Bodwin and Yennie for a calculation 
of the elastic terms only, using Eq. ([3J at the photon-proton vertices and no 
dispersion theory [9] . After adding the same term to the polarizability corrections 
in Ai, one obtains the commonly quoted result for A\ |15|16|18] , Beyond the 
historical connection, if one is comfortable with the unsubtracted dispersion 
relation, the use of the dispersion theory gives a more secure result because it 
uses only the pole part of the photon-proton-proton vertex, so that the combined 
elastic and inelastic result does not depend on the general validity of whatever 
photon-proton-proton vertex one uses. 

5 Conclusion 

The evaluation of the polarizability contributions to hydrogen hyperfine struc- 
ture, A po \, based on latest proton structure function data is A po \ = 1.3 ± 
0.3 ppm |Tj5j . This is quite similar to the Faustov-Martyenko 2002 result, which 
we think is remarkable given the improvement in the data upon which it is based. 
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Most of the calculated A po \ comes from integration regions where the photon 
four- momentum squared is small, Q 2 < 1 GeV 2 . 

There is still a modest discrepancy between the hydrogen hfs calculation and 
experiment, on the order of 2 ppm. Optimistically, one can hope for a rapid 
reconciliation between data and calculation. It surely has not been missed that 
using the Kelly form factor [13] value of the Zemach radius and the new Faustov- 
Gorbacheva-Martynenko value for A po i [31] give excellent agreement between 
theory and data. Nonetheless, one can argue that other choices are currently 
better. The integrals that give Az and Z\ po i emphasize the low-Q 2 region. The 
Sick form factors are the only modern ones that are tuned to fit best at low Q 2 , 
and the determination of Z\ po i in [19] is the only one that has explicitly used the 
lower Q 2 Jefferson Lab inelastic data. 

An interplay between the fields of atomic and nuclear or particle physics 
may be relevant to sorting out problem. For one example, the best values of the 
proton charge radius currently come from small corrections accurately measured 
in atomic Lamb shift [38] . Sick's value of the charge radius [13] , from the analysis 
of scattering data, is somewhat larger. The precision of the atomic measurement 
of the proton charge radius can increase markedly if the Lamb shift is measured in 
muonic hydrogen [39., which could happen in 2007, if the Paul Scherrer Institute 
accelerator schedule holds. In the present context, the charge radius is noticed 
by its effect on determinations of the Zemach radius. 

For ourselves, we look forward to a high accuracy resolution of the proton 
structure corrections to hydrogen hfs, and also to finishing a clear continuation 
of the present program by the evaluation of the muonic hydrogen ground state 
hfs. We have formulas with all lepton masses in place, and are currently waiting 
until the final EG1 data is released, which we think will be rather soon, before 
proceeding and publishing a numerical evaluation. 

My contributions to this subject have all been made in collaboration with 
Vahagn Nazaryan and Keith Griffioen. I thank them for the pleasure I have 
had working with them. In addition, we thank Jose Goity, Savely Karshenboim, 
Ingo Sick, Silvano Simula, and Marc Vanderhaeghen for helpful discussions and 
information. This work was supported by the National Science Foundation under 
grants PHY-0245056 and PHY-0555600 (C.E.C.); PHY-0400332 (V.N.); and by 
the Department of Energy under contract DE-FG02-96ER41003 (K.A.G.). 
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